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Abstract. Multi-flavour (Nf > 3) Chiral Perturbation Theory (xPT) may exhibit instabilities due to 
vacuum fluctuations of sea qq-pairs. Keeping the fluctuations small would require a very precise fine- 
tuning of the low-energy constants L 4 (a*) and L e (/i) to Lf lt (M p ) = -0.51 ■ 10~ 3 , £e rit (M p ) = -0.26- 10~ 3 . 
A small deviation from these critical values - like the one suggested by the phenomenology of OZI-rule 
violation in the scalar channel - is amplified by huge numerical factors inducing large effects of vacuum 
fluctuations. This would lead in particular to a strong 7V/-dependence of Chiral Symmetry Breaking (xSB) 
and a suppression of multi-flavour chiral order parameters. A simple resummation is shown to cure the 
instability of Nf > 3 X?T, but it modifies the standard expressions of some 0(p 2 ) and 0{p 4 ) low -energy 
parameters in terms of observables. On the other hand, for r — m s /m > 15, the two-flavour condensate 
is not suppressed, due to the contribution induced by massive vacuum ss-pairs. Thanks to the latter, the 
standard two-flavour xPT is protected from multi-flavour instabilities and could provide a well-defined 
expansion scheme in powers of non-strange quark masses. 



1 Introduction 

Understanding Chiral Symmetry Breaking (xSB) in low- 
energy QCD still deserves both phenomenological and the- 
oretical effort. First, there is a growing need to identify 
and to separate non-perturbative QCD effects from possi- 
ble manifestations of "New Physics" in experimental tests 
of the Standard Model (e.g. weak matrix elements, e'/e, 
(g — 2) M . . . ). Furthermore, the subject has its own theoret- 
ical interest. Vector-like gauge theories such as QCD for- 
mulated in a large Euclidean box allow a particularly at- 
tractive interpretation of spontaneous xSB in terms of the 
lowest modes of the Dirac operator averaged over all gluon 
configurations G). In QCD-like theories, some characteris- 
tic properties of the Dirac spectrum have been proven ^ 
|3| and possible consequences for chiral order parameters 
have been conjectured ji]. More generally, this approach 
to xSB suggests an ana ogy with disordered systems of 
higher dimensionality (d = 4) emphasising notions such 
as the average and the fluctuation of the density of small 
Dirac eigenvalues, as well as the transport properties (e.g. 
conductivity) [|| . Finally, the cornerstone of all this inves- 
tigation is Chiral Perturbation Theory (xPT) |||7]] which 
provides a systematic link between theoretical character- 
istics of xSB (described by order and fluctuation parame- 
ters) on one hand and observable properties of Goldstone 



bosons (masses, decay constants, scattering amplitudes, 
decay form-factors. . . ) on the other hand. 

During the last few years both theoretical and exper- 
imental progress was achieved along these lines. It has 
been suggested that order parameters of x^B, in partic- 
ular the quark-antiquark condensate (qq), could strongly 
depend on the number Nf of light flavours MM ■ As Nf 
increases, (qq) as well as the Goldstone boson coupling 
F„ are gradually suppressed, due to the paramagnetic be- 
haviour of Dirac eigenvalues and to increasing fluctuations 
of the density of states HQ]. This effect is induced by 
light-quark loops and it cannot be detected in quenched 
lattice simulations. Actually, there are two kinds of param- 
agnetic effects generated by loops of sea-quarks which are 
both of the same origin []: the massless loops suppress chi- 
ral order parameters whereas the massive sea quark pairs 
enhance them, as long as their mass is of order Aqc d or 
smaller. In Nature, this last remark merely concerns the 
strange quark, whose mass is slightly below Aq C d- The 
abundance of strange quark-antiquark pairs in the vac- 
uum can thus lead to a different behaviour of two-flavour 
(m s ~ Aqcd) and three-flavour (m s = 0) chiral dynam- 
ics. Such difference would be characterised by a nonnegli- 
gible vacuum correlation between strange and non-strange 
quark pairs, implying in turn large 1/N C corrections and 
violation of the OZI-rule in the scalar channel. The latter 
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is actually observed |L0j and a strong variation of xSB 
between Nf — 2 and Nf = 3 has been indeed reported 
on the basis of sum- rule studies @>ll>@ 

using as input 
available information about the scalar sector ++ . 

Such a possibility should now be considered in the 
light of the new experimental information on low-energy 
7T7T scattering which has been recently published Jl5[] and 

The outcome of these analyses 



analysed (T6|,|17|, |1S|, g_9j,|2_ 

shows that, in the presence of massive ss pairs in the vac- 
uum, the two-flavour condensate (uu) is large and domi- 
nates the SU (2) x SU (2) symmetry-breaking effects jl9[ 
|20|] . Accordingly, the standard two-flavour %PT expansion 
in powers of m u and m c i Q should be expected to converge 
rather well. On the other hand, large vacuum fluctuations 
of qg-pairs would result into a large difference between 
Nf = 2 and Nf — 3 condensates, destabilising the three- 
flavour expansion. Indeed, a detailed SU(3) x SU (3) analy- 
sis of Goldstone boson masses and decay constants within 
the standard two-loop x?T |2lJ has revealed an anoma- 
lously large 0(p 6 ) contribution to M%, depending on a fine 
tuning of the LEC's L±(n) and Lq(iS) - which precisely re- 
flect vacuum fluctuations. The purpose of this paper is to 
show that the difference in the chiral behaviour of two- 
flavour and multi-flavour QCD described above could be 
naturally explained in terms of the interplay between vac- 
uum fluctuations (of small Dirac eigenvalues) and chiral 
order (described by order parameters such as (qq)) ■ 

We start by considering Ward identities and low-energy 
theorems for the two-point functions (D a D b ) and (V a V — 
A a A b ), where V a , A a and D a are (charged) vector cur- 
rents, axial currents and the divergences of the latter, re- 
spectively. We write these identities in a form reminiscent 
of the xPT expansion of MpFp and Fp (P = n, K, rf), in- 
cluding explicitly the leading and next-to-leading orders 
in powers of quark masses |7j] and collecting (not neglect- 
ing) all remaining orders into well-defined "remainders" . 
We refer to Ward Identities written in this way as "mass 
and decay constant identities" . 

We then show that there exist two exact nonlinear 
relations between the order parameters (qq) and Fp in 
the chiral limit SU(N f ) x SU(N f ) with N f > 3 and two 
"fluctuation parameters" which are defined in terms of 
the standard LEC's Lq(^l) and L^/j,): in these relations, 
all the effects of higher x?T orders are absorbed into a 
finite multiplicative renormalisation of order and fluctua- 
tion parameters. The deviation from 1 of the correspond- 
ing renormalisation constants (rescaling factors) remains 
under control to the extent that the NNLO remainders in 
mass and decay constant identities are small. The x?T 
series is reproduced in the limit of small fluctuation pa- 
rameters, implying a very precise fine tuning of Lq((j.) and 
Li4,(p,). Otherwise (and in particular for large fluctuation 
parameters), the multi-flavour chiral condensate is sup- 
pressed, and the standard xPT interpretation of mass and 
decay constant identities breaks down. This need not af- 
fect the overall convergence of the xPT series. We actually 
expect that for the physical value of m s the multi-flavour 
XPT still makes sense globally. In this case the instability 
caused by large fluctuations and by the suppression of the 



quark condensate can be cured by a simple "resummation" 
which amounts to replacing the perturbative solution of 
the non-linear relations between order and fluctuation pa- 
rameters by their exact algebraic solution. This merely 
modifies the standard way of expressing the parameters 
of the effective Lagrangian in terms of observable quanti- 
ties §. 

The multi-flavour mass and decay constant identities 
can be further used to define the two-flavour order param- 
eters, by taking the limit m u , — > but keeping m s fixed 
at its physical value. In this way we show that the two- 
flavour condensate and decay constant are not affected 
by the large fluctuations which suppress the multi-flavour 
condensate. In addition, this allows one to discuss the con- 
nection with two-flavour xPT an d the new in: scattering 
data. 

The plan of this article is the following. In Sec. 0, we 
discuss the impact of loops of massive sea-quarks (typi- 
cally, ss-pairs) on the pattern of x^B. We introduce in 
Sec. I the mass and decay constant identities for a generic 
Nf > 3, and show that they lead to an exact system of 
relations between order and fluctuation parameters, pre- 
sented in Sec. || The general properties of this system are 
then discussed in the light of the positivity and (conjec- 
tured) paramagnetic inequalities that the order and fluc- 
tuation parameters have to obey (Sec. |J). Sec. ^is devoted 
to the study (in the plane of fluctuation parameters) of the 
critical line where the symmetry is restored, i.e. both the 
condensate and the decay constant vanish. Sec. fj] briefly 
summarises properties of the large- N c limit in which fluc- 
tuations are suppressed. Then in Sec. || we discuss the 
opposite limit of large fluctuations. A possible realisation 
of the latter can be interpreted as a limit of large Nf. It 
is shown that the limit of large fluctuations is in principle 
different from the symmetry restoration limit: despite the 
continuous vanishing of the multi-flavour quark conden- 
sate, the decay constant stays non-zero. Sec. ^| deals with 
the SU(2) x SU(2) chiral limit: we check once more that 
in the large-fluctuation limit the two-flavour condensate 
remains non-zero and the two-flavour Gell-Mann-Oakes- 
Rcnncr relation is approximately obeyed, and we discuss 
briefly the recent results on 7T7t scattering. The conclusion 
and a few appendices close the paper. 



2 Role of the mass of sea quarks 

When defining a chiral limit or a chiral order parameter, 
it should be stated which fermions are taken massless and 
which quarks in the sea are left massive. The simplest sit- 
uation is the one with just N massless fermions ip\ . . .tptf, 
(i.e. mi = mi = . . . tojv = m — > 0) and no other mas- 
sive fermions left. The condensate of this purely massless 
theory is defined as: 



S tj a(N) 



lim (ipitpj 

m— >0 



■N ■ 



(1) 



The expression of the condensate in terms of eigenvalues 
A„ of the Euclidean Dirac operator 7 M D AI [G] defined in a 
box L x L x L x L with periodic boundary conditions (up 
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to a gauge transformation) is well known as the Banks- 
Casher formula M. Formally, it can be written as: 



a(N) 



lim — <C 
V 



> 



where lim means taking V — > oo first and m — + after- 
wards. The average <C 3> over Euclidean gluon configura- 
tions involves the iV-th power of the fermion determinant 
A N (m,G). Since for m — > 0, V — > oo only the smallest 
eigenvalues contribute in Eq. (||) , it is conceivable that the 
main iV-dependence merely arises from the infrared part 
of the determinant A(m,G) — AirAjjv, defined as: 



A IR (m, G) = m 



n 

n<K 



m 2 + \l 



(3) 



where v is the winding number of the gluon configuration 
G and K(A, G) is an integer corresponding to a cutoff 
A such that \k — A. Positive eigenvalues are ranked in 
ascending order Ai < A2 < . . . < \k = A. The numbers 
uj n are defined by the Vafa-Witten bound || for the Dirac 
eigenvalues: 

/ ji \ l/d 

K<C (-) EE u> n . (4) 

The existence of such a uniform bound independent of 
gauge field configurations is a specific property of QCD- 
like gauge theories. It implies (Ajh) n+1 < (Air) n for 
any finite cutoff K and for m < fj,. In the chiral limit 
V — ► 00,777 — * one can then expect the paramagnetic 
inequality: 

<j(N + 1) < a(N) . (5) 

In the real world, the situation is slightly more compli- 
cated due to the role of massive virtual quark pairs which 
may be present in the vacuum. Notice that massive and 
massless pairs have different chiral transformation prop- 
erties and do not affect the chiral structure of the vacuum 
in the same way. In QCD one deals with a hierarchy of 
quark masses: 



m u < ITid -C m s <C ?77 c < 777/, <C 777 t . 



(6) 



Some of them (u, d, s) can be considered as light com- 
pared to the scale Ah ~ 1 GeV at which the masses of the 
first bound states non-protected by chiral symmetry oc- 
cur. Ah ~ ^wFir is the reference scale in x?T expansions 
in powers of p/ Ah- A different question can be asked in 
connection with the structure of the vacuum. Some quarks 
(c, b, t) can be considered as too heavy to form abundant 
vacuum pairs. In this case, the characteristic scale is not 
Ah but the lower scale Aqcd- The reason is that we are 
not interested in the production of massive hadrons but 
in the creation of massive virtual qq pairs. The latter will 
be most probable if the quark mass is of order Aqc d or 
slightly lower (if m — ► 0, the chiral properties of the cor- 
responding pairs tend smoothly to the massless case). 

This reasoning already singles out the strange quark 
among all six quarks we know. Its mass m s ~ 160 MeV 
(at n = 1 GeV) [0 is sufficiently low compared to Ah 



to legitimate a SU(3) x SU(3) chiral expansion. On the 
other hand, m s is sufficiently close to Aqcd to expect a 
significant presence of massive ss-pairs in the vacuum. 

We define the SU{N f ) x SU(N f ) chiral limit in QCD 
by taking the first Nf quarks as massless and keeping the 
remaining masses at their physical value. In practice, one 
can consider such a limit for Nt = 2 or for Nt = 3. The 
corresponding order parameters will be functions of the 
remaining non-zero masses. For instance, the two-flavour 
condensate is defined as: 



S(2,m s ,...) 



lim 

L,md- 



(7) 



and it is a function of m s as well as of the heavy-quark 
masses denoted in Eq. (0) by the ellipsis. u(x) stands for 
the lightest (u) quark field, and it can equivalently be 
replaced in Eq. (Q) by the d-quark, but not by the s-quark 
field. The three-flavour condensate is then defined as: 



2(3, 



.) = lim U(2, m s 

m s —>0 



(8) 



One expects that the effect of heavy-quark masses 777 c , mb, 
m t on 2(2) and ^(3) remains small and could be eventu- 
ally estimated. For simplicity, we shall neglect all effects of 
heavy-quark masses in the sequel. In this approximation, 
2(3) coincides with cr(3) defined in Eq. ([[]). It is a clean 
probe of the chiral structure of the vacuum of QCD with 
nothing but three massless quarks: once one sets ?77 s = 0, 
there is no more massive quark left which would be suffi- 
ciently light to pollute the vacuum \ f2)^. 

The situation is rather different in the two-flavour chi- 
ral limit m u = = m — > 0, keeping the strange-quark 
mass 777,5 at its physical value. Since 777 s is not very large 
compared to Aqcd and vacuum is polluted by massive 
ss pairs, it is difficult to relate 2(2) to the genuine con- 
densate a(2) characteristic of the theory with nothing but 
two massless quarks. This situation occurs as long as m s 
remains of order Aqcd- One can gain more insight into 
the ?77 s -dependence of ^(2) from the formula in Euclidean 
space: 

(uu exp[— 777 s J ss])s 



2(2) 



(9) 



(exp[-m s J ss}) 3 

where on the right hand side the expectation value is taken 
with respect to the vacuum of the SU(3) x SU(3) invari- 
ant theory. It is seen that in the absence of correlations 
between the strange and non-strange qq pairs, Eq. (Q) im- 
plies 2(2) = 2(3), in agreement with large-iV c expecta- 
tions. For rri s <C Ah, Eq. (||) can be rewritten as: 



2(2) = 2(3) + m s Z{ 



where Zf denotes the connected correlator of ss- 
pairs: 



lim f 



Z{ = lim / dx (ss(x)uu(0)) c 



(10) 
and uu- 



(11) 



Notice that in the limit m u ^ — * 0, uu can be equally re- 
placed by (uu + dd)/2. The correlator ([fl]) measures the 
violation of the OZI rule in the isoscalar scalar (i.e. vac- 
uum) channel and it can be estimated using the experi- 
mental information now available in this channel [hi. 12, 
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It is related to the standard 0(p 4 ) LEC L e (p) and it 
turns out to be larger than what is expected on the basis 
of large- N c considerations @ . 

It is useful to express E(2), E(3) and the correlator 
((ss)(uu)) in terms of the eigenvalues A„ of the Euclidean 
Dirac operator. Neglecting heavy quarks, both E(3) and 
E(2) concern the theory with the same total number of 
fcrmions: in E(3), one sets m u — — m s — to — > 0, 
whereas in E(2) m u = m^ = m — ► but m s ~ Aqcd 
is held fixed. Hence, the corresponding Banks-Casher for- 
mula becomes |BJ: 



lim-«£ 



(12) 



where the only difference lies in the determinant inserted 
in the average over gluon configurations: A 3 (m,G) for 
E(3) and A 2 {m,G)A{m s ,G) in the case of E(2). Com- 
paring the corresponding infrared parts Eq. (^) which are 
expected to dominate in Eq. flLq) , one has: 



m 



3\v\ 



n 

n<K 



<m a Mmjr'n (m ^5 ) , 2( ^, + A " ) . <™> 



n<_ff 



(M 2 + 



as long as to — ► and to s ~ Aqcd- This suggests that the 
paramagnetic inequality (^) holds even in the presence of 
massive strange quarks in the sea: 



E{2) > E(3) . 



(14) 



The conjectured inequality Eq. ( |14| ) will play an important 
role in the sequel. 

It remains to express the correlator Eq. ( |ll| ) in terms 
of Dirac eigenvalues; see e.g. Ref. M. One has 



Z{ = lim 



nk 



to 2 + A 2 to 2 . + Xl 



»c 



(15) 



For small to and m s , only small Dirac eigenvalues con- 
tribute and the expression ( )l5| ) measures their correla- 
tions. For m s r~j m, Zf describes the fluctuations of the 
density of states p(X) — J2 n ^(^ — ^™) near A ^ 0. The 
positivity of Eq. dig) is in agreement with the paramag- 
netic inequality Eq. (|1J). In this paper, we investigate the 
possibility that multi-flavour QCD vacuum behaves as a 
strongly correlated fermion-anti-fermion system charac- 
terised by large fluctuations Zf of the density of small 
Dirac eigenvalues. 

Let us mention without proof that the present discus- 
sion can be easily extended to the order parameter: 



F\N f ) 



lim 



F 



(16) 



3 SU(Nf) x SU(N f ) chiral symmetry for 

N f >3 

We consider QCD with Nt = n + 2 light flavours which 
is very like our actual three flavour QCD equipped with 
n copies of the strange quark degenerate in mass: u and d 
quarks with ultralight degenerate masses m u = to<j = to 
and n "strange" copies si...s n with a common mass m s 3> 
to but still light compared to the QCD scales. 

In the chiral limit to = m s = 0, the SU(2+n) x SU(2+ 
n) chiral symmetry is assumed to be spontaneously broken 
down to the diagonal subgroup SUy in + 2). This symme- 
try is explicitly broken by the mass difference m s — to to 
SU{2) x SU (n), which is the exact symmetry of our prob- 
lem. We have (n + 2) 2 — 1 = n 2 + 4n + 3 pseudo-Goldstone 
bosons: 3 pions of mass M % , An kaons (usi,dsi for i = 
1 . . . n, plus conjugates) of common mass Mr-, the 77-meson 
with structure A ?) = (l + 2/n)~ 1/2 diag[l, 1; -2/n . . -2/n] 
and mass M v , and finally the extra n 2 — l SiSj states whose 
mass will be denoted as Mx- 

This framework will be used to discuss the response 
of order parameters of the chiral symmetry SUz,(Nf) x 
SUji(Nf) to large vacuum fluctuations both for Nf > 3 
and for Nf — 2. We will see that these two cases behave 
rather differently. It will prove useful to keep a generic 
value for Nf, although our main interest will be on the 
case Nf = 3 and all our conclusions apply to this case. 



3.1 Multi-flavour mass and decay constant identities 

For all n we are using the same 0(p 4 ) effective Lagrangian 
as for the n = 1 case of Ref. (except for one additional 
invariant (V^U^VUV ^t/^V V U) , irrelevant for our pur- 
pose). The LEC's have an a priori unknown dependence 
on n, therefore we keep the same notation for all LEC's 
with an extra index n. This is justified to the extent that 
the symmetry properties are as described above. 

The chiral expansion of Goldstone boson masses is very 
close to the case n = 1 discussed in Ref. Q] : 

F 2 M 2 = 2mE(n + 2) + 2m(nm s + 2m)Z s n + Am 2 A n 

+Am 2 B 2 n L + F 2 M 2 d wn , (17) 

F%M\ = (to s + m)E(n + 2) + ( m s + m)(nm s + 2m) Z^ 
+ (m s + m) 2 A n + m(m s + m)B 0n L 
+m s (m s + m)B 2 n L' n + F 2 K M 2 K d Kn , (18) 

F X M X = 2m s E(n + 2) + 2m s (nm s + 2m)Z s n + 4m 2 s A n 
+4m 2 B 2 0n L' n + F 2 M 2 x d Xn , (19) 

and likewise for F 2 M 2 (see App. ^J). The connection with 
the standard LEC's of the Nf > 3 effective Lagrangian is: 



defined as before by taking the first Nf quarks massless 
and leaving the remaining quarks at their physical masses. 
The dependence of these order parameters on Nf and on 
the sea-quark masses is qualitatively similar to that of the 
quark condensates. 



32B 2 n iLftn) 



5127T 2 



, M K 2 
log -I 

p? (n + 2) 2 



M 2 
log^- 
M 



(20) 
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A n = 16B 2 J 



5127T 2 

We have the ratio: 



, Ml 



M 2 

-^log^- 

2 /i z 



Bon 



S(n + 2) 
F 2 {n + 2] 



(21) 



(22) 



and E(n + 2) and F(n + 2) denote the condensate and 
the pseudoscalar decay constant for all n + 2 massless 
quarks. The remainders FpMpd 2 P collect all higher-order 
terms, starting at the next-to-next-to-leading order O(m^) 
(NNLO) [hence d n — 0(m 2 )]. The (n-independent) com- 
bination of chiral logarithms L has the same meaning as 
in Rcf. Ill: 



3.2 Low-energy constants 

It is a simple exercise of algebra to combine Eqs. (p7[)-(p^|), 
and Eqs. (|25|)-(|2^), to get the expression of the LEC's: 



2m 



[E(n + 2) + (nr + 2)mZ* 



1 - e(r) - d„ - 



4m 2 B 2 n 



F T 2 M T 2 r-1 



(L-L^)(30) 



im 2 A n 
JfM 2 



e(r) + d' n 



4m 2 B 2 n 1 



— (i-ri^), (31) 



ith 



e(r) 



> z 

H — 1 



F K M K 
F^U 



1 



(32) 



d n and ci^ are linear combinations of the remainders d„ 
and dKn- 



L 



32tt 2 



o, M 2 K M 2 

3 log — -£r + log — 

M 2 m 2 



A' 



(23) 



whereas L' contains all chiral logarithms involving the un- 
physical mass Mx- 



L' = 



— |log^ 
16^ 2 n } 8 M 2 K 



(n + l)lo{ 



M\ 



(24) 



which vanishes for n — 1 as it should. Similar expressions 
are derived for the decay constants: 



F 2 = F 2 (n + 2) + 2(nm s + 2to)£„ + 2m£ n 
2mB 0n 



(25) 



+ 



16tt 2 



M 2 

l0g 4 



A/ 2 



F 2 e 



F| = F 2 (n + 2) + 2(nm s + 2m)£ n + (m s + m)£ n (26) 



+mB 0n L + m a B 0n L' n + F%e Kn , 
F 2 (n + 2) + 2(nm s + 2m)£„ + 2m s £„ (27) 

, 2m s B 

- 1 - ""*."T iK 



while F v is given in App. |A[ The scale-invariant constants 
£ n and are related to the LEC's L4 and L5 as follows: 




256tt : 



M 2 
Tog— 



(28) 



nlog^ + 21og^ 



(29) 



The remainders Fpep n collect the higher-order terms, start- 
ing at NNLO [0(m 2 )}. 



r+L / , , 2 
-d„ n - e(r) + 



dKn 1 d n — dn d-ffn ■ 

(33) 



r — 1 \ r — 1 

A similar work can be performed for the decay constants: 

2mB 0n 



- rj(r)+e' 



F 2 



n (r~l)F 2 



L + tL> 



1 / Ml M 2 
"ft? l0g 4 + 21 ° g M| 



(34) 



. 2m£ 

(nr + 2)-pj- = 1 - - e, 



2mB 0n 




F 2 {n + 2) 
F 2 

L + rL' 



M 2 



2 log 



K 



M 
Ml 



K 



,(35) 



ith: 



r + 1 . . 2 

e n = re™ - r/(r) + 



r- 1 



r- 1 



(36) 

&Kn = e' n + e nn (37) 



We recall that e(r) and r/(r) are suppressed at large values 
of r (> 15). We expect then d' n <§C d n ~ and <C 

One easily checks that all the formulae displayed in 
this paragraph reduce to the mass and decay constant 
identities obtained in Ref. |14| in the case n = 1. 



3.3 ,577(2) x 577(2) order parameters 

In the SU(n + 2) x SU(n + 2) chiral limit (m = m s = 0), 
we define the dimensionless order parameters: 



X(n + 2) 



2mS(n + 2) 
F 2 M 2 
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Y(n + 2) 



Z(n + 2) 



2mB, 



Oh 



Ml 



(38) 



F 2 (n + 2) 
F 2 



The Gell-Mann-Oakes-Renner ratio X(n + 2) measures 
the quark condensate in physical units, while Z(n + 2) 
does the same for the decay constant. We have X(n + 2) = 
Y(n + 2)Z(n + 2). 

We consider now the SU (2) x SU(2) chiral limit where 
only m vanishes (and the n copies of the strange quark re- 
main massive), in order to investigate the effect of massive 
sea quarks on two- flavour chiral dynamics. We define the 
quark condensate and the decay constant in this limit as: 



E{2) = lim 



F1K 



F 2 {2) = lim F 2 



m— >0 2l7l ' v ' m— >Q 

and the dimensionless order parameters: 
2mS{2) 



(39) 



X(2) 
Z{2) = 
Y{2) = 



F 2 M 2 ' 
F 2 {2) 

7T 

X{2) 2mB{2) 



(40) 



Z(2) 



Ml 



From Eqs. (|17]) and (p0|), we derive the expression for 
the two-flavour condensate: 



X(2)[l-4n] = 



T! T 



nr + 2 
2 



1 - e(r) - d n - [Y{n + 2)] 2 /« 
-X(n + 2), (41) 



where a quantity topped with a bar is considered in the 
limit m — ► and: 



fit 

K 

AZi = 



Ml 

n 

Zi H 



16tt 2 



B 0n AZ n ) 
A" 



M 



2 M 2 

lOE - 

(n + 2) 2 S M 2 



(42) 



(43) 



In a similar fashion, the two-flavour decay constant can 
be read from Eqs. (|25|) and (^5|): 



Z(2)[l 

with the quantities 

M 2 ( ■ 



nr 



nr + 2 
2 



+ - 



nr 



[1 - n{r) - 
-Z(n + 2). 



Y{n + 2)g v 



(44) 



-{nr + 2)AC 



(45) 



/„ and g n contain chiral logarithms of two different types 
and signs: the first involve the masses of the Goldstone 
bosons M^^k,t),x and are positive, while the latter are 
negative and combine ratios of masses of the same meson, 
considered in the massive theory and in the SU(2) x SU (2) 
chiral limit (m — > 0). We can estimate /„ and g n by iter- 
ating the previous mass and decay constant identities as 
explained in Apps. ^ and |b[ 



4 Connection between fluctuation and order 
parameters 

We want now to investigate the connection between chiral 
order parameters and vacuum fluctuations. For a generic 
n > 0, the mass and decay constant identities yield a very 
simple expression of the two order parameters of funda- 
mental interest, namely the condensate S{n + 2) and the 
decay constant F(n + 2), in terms of the two OZI-rule vi- 
olating LEC's L 4 and L 6 . These two large- N c suppressed 
constants directly reflect vacuum fluctuations of qq pairs 
and they provide a convenient framework to discuss com- 
pletely and transparently how these fluctuations affect or- 
der parameters. 

The first step consists in eliminating the LEC's A n and 
£ n from the mass and d ecay constant identities for the 
pion and the kaon Eqs. and leading 

to Eqs. ( |30| ) and (|35|). Then we can reexpress the (OZI- 
rule violating) constants £ n and Z^ in terms of L4 and 
Lq, and obtain the two desired relations between the order 
parameters X, Y, Z and the fluctuation parameters L4 and 

X(n + 2) = 1 - c (r) - d n - [Y(n + 2)} 2 p n /4 , (47) 
Z(n + 2) = 1 - n(r) - e n - Y(n + 2)A„/4 . (48) 

The constants Lq and L4 enter the discussion through 
the combinations: 

Pn = 64 ^f>r + 2 ) AL 6 > A„ = 32^f (nr + 2)AL\ , 

(49) 

where the scale-independent differences AL™ = L"(/i) — 
L"' CTlt (n) involve the critical values of the LEC's defined 
as: 



r n,crit/ \ 1 , M K 



256tt 2 a 2 8(r-l)(nr + 2) 



[L-L'n] 



1 r+1 1 M V 

1o Zt£> ( 5 °) 



256tt 2 (r- \){nr + 2) M 2 



6 (W = 



1 



M 



M 2 

2 1°§ ~2~ 



512tt 2 V b u 2 (n + 2) 2 b n 
[L-L'J. 



16(nr + 2)(r- 1) 



(51) 



32tt : 



log 



Ml 

Ml 



(46) 



Comparing the fluctuation parameters A n and p n ( |49| ) 
to 1 provides a quantitative measure of the effect of the 
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LEC's L4 and Lq on observable quantities. The effect dis- 
appears if La and Lq are fine-tuned to their critical values 
( p(\ ) and (|5l|), which for n = 1, r = 25, fj, = M p become 
Lf lt (M p ) = -0.51 • 10~ 3 and L c 6 Tit (M p ) = -0.26 • 10" 3 . 

Notice that even a small deviation from the critical 
values is amplified by the large numerical coefficients in 
Eq. (^) . ft is possible and convenient to absorb the NNLO 
and higher-order contributions - represented in Eqs. (|47|)- 
( ff8| ) by the remainders d n and e n - into a multiplica- 
tive renormalisation of order and fluctuation parameters. 
Defining the rescaled parameters as: 



F 2 (n + 2) 



i^[l-r?(r)-e„]' 



Dn 



1 - r)(r) - e n 2mB 0r , 



l-e{r)-d n Ml 



(52) 



Wn — ^nkn> V n — Priori > ^n(^*) 



1 - e(r) - d n 
(1 - i](r) - e n ) 2 



(53) 



the two basic equations inferred from the mass and decay 
constant identities take now the concise form: 



1 

1 1 

4 y, 



(54) 
(55) 



In addition, the multi-flavour quark condensate [ex- 
pressed in GOR units, cf. Eq. Q3q)] can be rewritten in 
terms of y n and z n : 



X(n + 2) = [1 - e(r) - d n ]y n z n 



(56) 



At this place, a few remarks are in order: 

i) The above analysis holds for a generic n > 0, in- 
cluding the physical case n = 1 of three flavours. On the 
other hand, the case of two massless flavours (m = 0, m s 
nonzero) requires a separate discussion based on Eqs. ( ^l| ) 
and (HJ) above. 

ii) The relations (IJ) and (H|) between the (rescaled) 
order parameters y, z on one hand and the fluctuation pa- 
rameters u, v on the other hand represent exact identities 
which do not result from any approximation or expansion. 
The influence of higher chiral orders - NNLO and beyond 
- is entirely encoded in the rescaling factors (^) and ( |53"|) 
through the remainders d n and e„. The latter, defined in 
Eqs. ( [33] ) and (|37|), stem from 0(p e ) (and higher) contri- 
butions to the mass and decay constant identities. They 
are both of order d n ,e n = 0(m 2 ) and should remain small 
unless the whole xPT series blows up. 

in) We expect the rescaling factors in Eqs. ( p3^ ) and 
( |53| ) to be close to 1. Two circumstances could spoil this 
expectation. First, if the quark mass ratio r were small 
(typically r < 10), the quantity e(r) would get close to 
1. On the other hand, for r > 15 one has e(r) < 0.2 and 
r](r) < 0.07. Such higher values of r are preferred by the 
new low-energy 7T7T scattering data ^0|, discussed in Sec. ^|. 
Second, even if the remainders d n and e„ are small in the 
physical case n = 1 (say 10 % or less), their size could 
grow when n increases. We shall come back to the large-n 



behaviour of the theory in Sec. |^. In either case, we would 
not be allowed to treat perturbatively the rescaling factors 
as close to 1. 

iv) Let us consider a solution of the generic system (|54|) 
and ( |55| ) giving the order parameters y and z in terms of 
the fluctuation parameters u and v. We can then read 
from Eqs. |n]) and @ together with Eqs. @ and @ 
the following parameters of £ e g: 2mB 0n , F§ = F 2 (n + 2), 
L^(p.) and Lg,([i) Qas functions of r = m s /m, L^(fJ,), Lq(h) 



the NNLO remainders d^m dx-a and e nn , ex-n Q These ex- 
pressions are then used in the study of other observables 
(e.g. Goldstone Boson scattering amplitudes) to eliminate 
0(p 2 ) and 0(p 4 ) LEC's from the bare xPT formulae. In 
particular, this procedure should be used when the NNLO 
remainders are (iteratively) matched with two-loop xPT 
expressions. The procedure described above - eliminating 
constants of C e s in favour of observables - differs from 
the one used in standard xPT, unless the fluctuation pa- 
rameters u and v are small compared to 1. This might in 
particular affect the outcome of standard analyses beyond 
one loop |pl[ . 

v) If u n ,v n <C 1, the standard multi- flavour x?T may 
be recovered in two steps. First, one constructs the per- 
turbative solution of the non-linear system (|Z|) and (IBT 
in powers of u n ,v n = 0(m q ): 



Vn = 1 + ^(«n - V„) + i(w„ - V n ) 2 + 



1 - \u r , 



1 

16 



u n (u n - v n ) + . . 



(57) 
(58) 



Next, one returns to the original variables A„, p n in 
Eq. @ and to 2mB 0n /M 2 F 2 (n + 2)/ F 2 expanding the 
rescaling factors in Eqs. <p2\ ) and ( |53|) around 1 in pow- 
ers of e(r) = 0{ni q ), rj{r) = 0(m q ) and of the remainders 
dn,en — 0(m 2 ). Matching the latter with the explicit two- 
loop contributions to Eqs. and @- (§(]), one re- 
produces the standard xPT expansion up to and including 
0(p 6 ) order. 

vi) In the physical situation of three massless flavours 
(n = 1), the fluctuation parameters need not be small 
enough to allow the power expansion of Eqs. (p7|)-(pq): 
using the estimates for L\{M p ) = 0.6 • 10~ 3 (central value) 
based on sum rules for the correlator {{uu + dd)ss) [ p] , |Ta 
[Tif and the recent determination of Ll(M n ) = 0.1 • lCf^ 3 
(central value) from ttK scattering data |2j| , one obtains 
the rough estimate u\ ~ 1.2, v\ ~ 3.4 (for r = 25). As 
pointed out in Ref. B, important vacuum fluctuations of 
q<7-pairs suppress the three-flavour condensate and desta- 
bilise the X-PT expansion. 

This phenomenon is a particular consequence of Eqs. 
(|54|) and (|55|). Multiplying the latter by y and using Eq. 



2 If we include the identities for r/, we can also express the 
constant L7 - see App. 

3 Alternatively, we can express all 0(p 4 ) symmetry breaking 
LEC's as functions of r and the two order parameters X(n + 2) 
and Z(n + 2), as discussed in Ref. jl4| for n = 1. 
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one obtains the relation: 



X(n + 2) = 2 



1 - e(r) 



1 + y/1 + V n /Z 2 n 



(59) 



which is identical (in the physical case n = 1) to Eq. (9) 
of Ref. || . 

To cope with possibly large values of the fluctuation 
parameters u and v, it may very well be sufficient to re- 
place the perturbative solution (|57|)-(|5^) of the system 
(|5^)-(|55|) by its exact algebraic solution, and to keep the 
perturbative expansion of the rescaling factors in Eqs. ( |52|) 
and (|53|). Let us emphasise that the divergence of the 
power series (|7])-([38|) is a question logically disconnected 
from the convergence of the expansion of rescaling factors 
around 1: the former is related to a possible suppression of 
multi-flavour condensate and of the corresponding leading 
order of xPT, whereas the latter is more a question of a 
global convergence of the x?T series starting at NNLO 
order. 

In the remaining sections of this paper we concentrate 
on the non-perturbative analysis of the system (p3)-(pq) 
and its consequences for the breaking of both SU{n + 2)x 
SU (n+2) and SU (2) x SU{2) chiral symmetries in various 
limits. 



5 Positivity and paramagnetic constraints 

In this section we investigate the system of equations (^)- 
( EE ) in the light of the positivity of the order parameters 
X(n + 2), Y(n + 2) and Z(n + 2) and of the conjectured 
paramagnetic inequalities they have to satisfy 0] . We are 
mainly interested in the domain allowed by these con- 
straints in the plane of the fluctuation parameters u and 
v. For simplicity, we omit further reference to n and r, 
since the latter enter the system (|5^)-(|5^) only via rescal- 
ing factors. The system admits two solutions for the ratio 
of the order parameter y (related to Bo n ): 



i + VTT" 



y- 



l - vT+u 



(60) 



These two solutions depend on the difference of the two 
rescaled fluctuation parameters u (function of the LEC 
L^) and v (function of L§). Notice that v — u is related to 
the particular combination 2L§ — L±. The square root in 
Eq. ( p0[ ) signals a non-perturbative resummation of vac- 
uum fluctuations as discussed in the previous section. For 
\v — u\ <C 1, the Taylor expansion of y + reproduces the 
xPT series, whereas ?/_ is of truly non-perturbative na- 
ture. 

The two branches y+ and y_ can be considered as two 
different sheets of the two- valued function y(v,u). These 
sheets are tangent along the line v — u = —1 that coincides 
with the boundary of the definition domain v > u — 1. 
Along this boundary, y has the constant value 2. The two 
branches y+ and y- are drawn in Fig. ^ as functions of 
v — u. When v — u increases, y+ decreases and vanishes 



asymptotically at infinity, whereas y_ increases and tends 
to infinity for v — u — > — . 

The system (|5^)-(|55|) yields the ratio y and the (rescaled) 
decay constant z, both of which should be positive for 
an acceptable solution: the vacuum stability requires x — 
yz > and z is related to F 2 (n + 2)/F 2 . As far as y is con- 
cerned, is positive in the whole half-plane v — u > —1 
(where < y < 2), and y_ is positive inside the strip 
— 1 < v — u < (where y > 2). The positivity of z, i.e. the 
condition z = 1 — uy/4 > 0, yields additional constraints 
in the (u,v) plane, especially for u > 0. The critical line 
z = is the parabola v = m 2 /4, along which y = 4/u. The 
condition z > is trivially satisfied for negative u, but 
for u > it leads on both sheets to (different) non-trivial 
bounds. 

As a result the whole domain of positivity of both y 
and z is obtained. On the y + sheet, one must have v > it— 1 
for u < 2, and v > u 2 /4 for u > 2. On the y- sheet, 
the positivity domain amounts to the part of the strip 
u > v > u — 1 situated below the parabola v — u 2 /4. 
These two domains are represented in Fig. |^. 
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6- 



y+ 
y- 



-0.5 



0.5 

v-u 



1.5 



Fig. 1. The two branches for y as functions of v — u. The solid 
(black) branch is y+, while the dashed (red) branch is 



We now arrive at the constraints imposed by the para- 
magnetic inequalities between order parameters discussed 
in Ref. j|. It was suggested that the chiral order param- 
eters that are dominated by the lowest eigenvalues of the 
Euclidean Dirac operator are particularly sensitive to a 
paramagnetic suppression arising from the infrared part of 
the fermion determinant - i.e. from the light-quark loops. 
This applies in particular to the decay constant F 2 and to 
the quark condensate S: the more flavours from the total 
of Nf = n + 2 become massless, the more suppressed the 
fundamental order parameters get. As a corollary, S(2) 
should be an increasing function of the strange quark mass 
m s , as long as the latter is comparable to the QCD scale. 
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Fig. 2. The two sheets for y in the (u,v) plane: on the left, 
the y+ sheet, on the right, the t/_ sheet. The hatched (green) 
regions are the allowed domains where y and z are both posi- 
tive. The value of y at each boundary is indicated. The critical 
line z = of symmetry restoration is denoted with a dashed 
(blue) line. 



This suggests the two paramagnetic constraints: 
X(2)>X(n + 2), Z{2)> Z(n + 2), 



(61) 



which can be reexpressed using the expressions of X(2) 
in Eq. @ and Z{2) in Eq. (||). They can be further 
simplified using the notation of the previous section, and 
yield for X and Z respectively: 



X : v > 



1 -D 



kf 



y 



Z :u> 



kg 



E 

y 



and 



1 - E 

' (62) 

where the NNLO remainders are rescaled: D n = d ml -(nr- 
2)/(nr), E n = e ln ■ (nr+2)/(nr), and / and g are positive 
combinations of chiral logarithms defined in Eqs. 
( filf ) and estimated in App. 

The paramagnetic inequalities Eq. ( |6l] ) yield there- 
fore lower bounds for the fluctuation parameters u and 
v. These bounds are only useful if the ratio of order pa- 
rameters y is large enough, for instance on the second 
sheet where y could grow. But we do not expect Eq. ( |62|) 
to have much relevance, say, on the first sheet for large u 
and v. Using the system of Eqs. (|54|)-(|55|), we may con- 
vert these lower bounds on the fluctuation parameters into 
upper bounds for y. The constraints on X and Z yield re- 
spectively: 



X 



y 



< 



Z:y< 



(l-D)z + y/(l-D) 
I- z{\- E) 

kg 



4k f 



(63) 



From Eq. (|63|) and the discussion of the previous sec- 
tion, we obtain the important result that y is necessarily 
non- vanishing and finite on the physical domain of the two 
sheets: 

< y < min ( -}=, ^- ) . (64) 



W k 9j 

Moreover, if y is large enough, Eq. (|6^) shows that u is 
positive, and thus z = 1 — uy/A < 1. The paramagnetic 
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0.5 



-0.5 



-1.5 





-1 
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0.5 
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5 2 



Fig. 3. Illustration of the paramagnetic inequalities on the two 
sheets for n = 1 and r = 25 (upper: first sheet y+, lower: second 
sheet y- ) . The boundaries of the two sheets are indicated with 
solid (green) lines, while the (blue) dashed line is the critical 
parabola z = 0. The (black) line with filled circles indicates the 
lower bound for v, derived from the paramagnetic inequality 
for X. The (red) line with crosses indicates the lower bound for 
u, derived from the paramagnetic inequality for Z. The NNLO 
remainders D and E are set to 0. 



inequalities in Eq. (|6l|) lead therefore to bounds for the 
two main SU (n + 2) x SU(n + 2) chiral order parameters 
X and Z. 

In the case n = 1 and r — 25,the estimation procedure 
detailed in App. [b| leads to Fig. || for u and v. We see that 
only the upper right part of the (it, v) plane survives. In 
particular, only a small fraction of the second sheet, far 
from the origin, remains available. Let us emphasise that 
the paramagnetic inequalities yield such constraints only 
if the two combinations of chiral logarithms / and g are 
positive. 

It is worth noting that the existence of minimal values 
of fluctuation parameters u and v does not contradict x?T 
which requires u, v — > in the chiral limit. The minimal 
values ( |6"2| ) are indeed proportional to M% as the parame- 
ters u and v themselves, and therefore constrain only how 
quickly u and v vanish in the chiral limit. 

As a conclusion, the system (|5^)-(^5|) admits two solu- 
tions for the ratio of order parameters y: they can be con- 
sidered as the two sheets of a two- valued function y(v,u). 
In the (u, v) plane of fluctuation parameters, the domain 
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of definition of the two sheets is restricted by positivity 
constraints [vacuum stability F 2 (n + 2) > 0] and by para- 
magnetic inequalities [X(n + 2) < X(2), Z(n + 2) < Z{2)\. 
As a result, the first sheet is limited to positive (and pos- 
sibly very large) values of the fluctuation parameters u, v, 
whereas only a tiny region of the second sheet, far from 
the origin but below u = 2, v = 1, is allowed. 



6 The case of symmetry restoration 

We investigate now the vicinity of the critical line where 
the SU(n + 2) x SU(n + 2) chiral symmetry is restored. 
Since we exploit results from the effective Lagrangian de- 
scribing the spontaneously broken phase in terms of Gold- 
stone bosons, we can comment only on the approach to 
the chirally symmetric phase from the broken one (and 
not on specific properties of the former phase, where the 
SU(n + 2) x SU(n + 2) chiral symmetry is restored). 

The symmetry restoration is equivalent to the van- 
ishing of the decay constant F„ in the chiral limit, viz. 
Z(n + 2) = 0. fn the (it, v) plane, the condition z = cor- 
responds to the critical line v — u 2 /4. The limit z — > + 
can be reached on both sheets: on the y+ sheet the critical 
line is approached from above and u > 2, while on the y_ 
sheet, the parabola v = u 2 /4 must be approached from 
below and < u < 2. On both sheets, the ratio of order 
parameters has the value y — 4/u on the critical line. 

We see that the vicinity of the origin in the (it, v) plane 
on the second sheet plays a special role in the discussion: 
y diverges there. The paramagnetic inequalities Eq. ( |64[ ) 
however prevent us from reaching the vicinity of it = v = 
on the second sheet: y is bounded and cannot become 
arbitrarily large. The ratio of order parameters Y(n + 2) 
remains thus finite when the critical line Z(n + 2) = 
is approached, so that the quark condensate X(n + 2) = 
Y(n + 2)Z(n + 2) vanishes. This was expected since all 
SU(n + 2) x SU(n + 2) order parameters must vanish as 
chiral symmetry is restored. 

It is also interesting to study the order parameters de- 
fined in the SU(2) x SU{2) chiral limit m = 0. Eqs. © 
and ( p[ ) yield the two-flavour decay constant and quark 
condensate: 



Z(2)[l 



X(2)[l-d nn ] = 



nr 


nr + 


2 


4 




nr 




nr + 


2 


4 




nr 




nr + 


2 



[1 - r)(r) - e 

y n kn(r)g n 4 

[1 - T)(r) - e 
4fe„(r)g„ 



nr 



(66) 



[1 - e(r) - d n ] 



x{l- yik n {r)f n + —y n z n \{Q7) 
nr 



16k n (r)f n 



(68) 



As long as we are not in the vicinity of the origin 
u = v = 0, we expect X(2) and Z(2) to remain close 
to 1 (for r > 20) along the critical line where the SU (n + 
2) x SU (n + 2) order parameters vanish. We can under- 
stand it by combining the definition of X(2) and Z(2), 
see Eq. ([40|), with the Ward identities for the masses and 
decay constants: 



X(2)[l-d m ] =X{n + 2) + nr 



2m 2 Z s n 
F 2 M 2 ' 



Z{2)[1 ~ e m ] = Z(n + 2) + nr 



2m^ 
F 2 



(69) 



(70) 



z „_>o nr 



nr + 2 



[1 - e(r) - d n ] 



For instance, the two-flavour quark condensate is the sum 
of its SU (n + 2) x SU(n + 2) counterpart and a LEC de- 
scribing the violation of the OZI rule in the scalar sector. 
We shall call the first term the "genuine" condensate - 
stemming directly from the breakdown of SU (n + 2) x 
SU(n + 2) chiral symmetry - and the latter the "induced" 
condensate - induced by the massive strange- like quark 
pairs present in the SU (2) x SU(2) vacuum [M. The same 
analysis applies to the decay constants Z(2) and Z(n + 2). 
We see now clearly how SU(2) x SU(2) chiral symme- 
try can remain broken while the SU (n + 2) critical line 
is reached. Even though there is no genuine contribution 
to the two-flavour order parameters, the vacuum is not 
invariant under SU(2) x SU(2) chiral rotations because of 
the symmetry-breaking transitions SiSi — > uu + dd which 
violate the OZI rule and are suppressed for N c — > oo. 

Let us now move along the critical line towards the 
origin u = v — by changing the value of the fluctuation 
parameters u and v (we must be on the second sheet to 
do so). Eqs. ( |66| ) and d68| ) indicate a suppression of X(2) 
and Z(2), leading finally to their vanishing: 

Z{2) > <-► u > 4kg, X{2) > <-> u > , (71) 

Which order parameter vanishes first depends on the rel- 
ative size of / and g: if g > y~fjk, Z(2) vanishes before 
X{2) does (and the other way round otherwise) . The point 
where at least one of the conditions Eq. ( [7l| ) is fulfilled 
marks the endpoint of the critical line : if we could pro- 
ceed further down the critical line, we would end up with 
an unphysical situation where X(2) or Z(2) is ne gati ve. 

One can check that the two conditions in Eq. ( 7l| ) are 
equivalent to the paramagnetic inequalities Eq. ([62 ) along 
the critical line z = Q - where y — 4/it = 2/^/v and 
SU(n + 2) x SU(n + 2) order parameters vanish. This 
was expected, since the paramagnetic inequalities along 
the critical line yield X(2) > X(n + 2) = and Z{2) > 
Z(n + 2) = 0, i.e. reduce to positivity constraints for the 
two-flavour order parameters. On the second sheet, the 
critical line ends thus at the edge of the domain allowed by 
the paramagnetic inequalities. In the physical case n = 1, 
and for r = 25, the right hand side of Fig. || shows that 
the inequality for v (i.e. for X) is saturated first and that 
the endpoint of the critical line corresponds to X(2) = 0, 
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Z(2) > 0. From App. g we see that this occurs in the 
physical case for any (large) value of r. 

In this section, we have investigated the critical line of 
SU (n + 2) x SU (n + 2) symmetry restoration. Along this 
critical line, both Z(n + 2) and X(n + 2) vanish, while 
Y(n + 2) = X(n + 2)/Z(n + 2) remains non-zero and fi- 
nite. We have studied the two-flavour order parameters 
X(2) and Z(2) as well. In the regions of the two sheets al- 
lowed by the paramagnetic inequalities, both are different 
from 0. There is only one exceptional point on the second 
sheet, where the critical line and the most stringent para- 
magnetic bound intersect. At this endpoint of the critical 
line, one of the two SU(2) x SU(2) order parameters van- 
ishes - X(2) in the physical case n = 1. 



7 No-fluctuation limit: N r 



OO 



We briefly sum up in this section the properties of another 
interesting region in the (u, v) plane where the effect of 
vacuum fluctuations is suppressed. This can be realised as 
the large- N c limit of the theory. Since 



u = 0(1/N C ) , v = 0(1/N C ) 



(72) 



we deal with the vicinity of the origin on the first sheet. 
The large- iV c limit forces the perturbative solution of the 
generic system Eqs. (|54|)-(|55|): 



y = l + 0(1/N C ) , z = l + 0{l/N c ) 



(73) 



The analogy between 1/N C expansion and xPT ceases 
here. The large- N c limit does not have much to say about 
the expansion of the rescaling factors in Eqs. ( |52| ) and ( |53"| ) 
- except perhaps by providing estimates of higher-order 
counterterms included in the NNLO remainders d n and 
e„. Combining the result ([73]) with Eqs. (jf^) and ( p3| ) one 
gets for large N c : 

X(n + 2)^1- e(r) - d n , Z(n + 2) -> 1 - v (r) - e„ , 

(74) 

From Eqs. and (f44|), we are now able to infer the two- 
flavour order parameters X{2) and Z(2) in the large-Af c 
limit. One of course expects to reproduce Eq. ((74|) - any 
dependence on the number of flavours should disappear at 
the leading order in l/N c . One first remarks that: 



/„ = 0(1/N C ) , g n = 0{l/N c ) . 



(75) 



so that for N c — > oo the right hand sides of Eqs. (H) 
and (|i|) both coincide with the one given in Eq._ (|74|). 
It remains to show that the NNLO remainders d vn = 
lim m ^o d^n and e„ = lim m ^o e^n are also suppressed 
for large N c . Indeed, the large-A c counting of connected 
QCD correlation functions yields: 

l im ^l og(j p2 M 2 )=0i Um log(F 2 )=(); 



Af c — >oo dm s 



Af c — >oo dm s 



(76) 

order by order in powers of quark masses, since the in- 
volved quantities receive contributions only from connected 



graphs with two (and more) fermion loops. This shows the 
suppression of d ml and e^ n and thus the desired result: 

X(n+2) = X(2)+0{l/N c ) , Z(n+2) = Z(2)+0(1/N c ) . 

(77) 

Both the condensate and the decay constant are indepen- 
dent of the number of massless flavours, and for not too 
small r both X and Z are close to 1. For completeness, 
we mention the leading large- N c behaviour of the LEC's 
Ls(/i) and Lg(fi), which can be read off from Eqs. (|31 
and d34|): 



L 5 = 8§3fo(»0 + ~ V(r) - e n ]k n (r) , (78) 



F, 



16M? 



z-[e(r)+d' n }[l-e(r)-d n }k n (r). (79) 



Notice that the scale dependence only shows up at the 
next-to-leading order. The OZI-rule suppressed constants 
1/4 and Lq remain 0(1) at large A^ c . They determine how 
quickly the fluctuation parameters vanish as iV c increases. 



8 The limit of large fluctuations 

We have seen that both multi-flavour (Nf > 3) xPT and 
I/A^ c expansion require small fluctuation parameters u 
and v. On the other hand, in the real world with three 
light flavours (n = 1), the fluctuation parameters Ui, t>i, 
as well as the difference _ui — v\ , are not small compared to 
I. Sum- rule studies fl]][l^ ,p^| suggest a significant though 
small deviation of the LEC's Ls(u) and L^fi) from their 
critical values Eqs. (fi7j|) and (|5l|). The large coefficients 
in Eq. ( ff9|) then amplify this deviation resulting into fluc- 
tuation parameters well above I. We will now study the 
effect of large fluctuations on the pattern of chiral sym- 
metry breaking. 



8.1 Many-flavour limit 

QCD does not contain an obvious parameter which could 
allow one to describe the limit of large fluctuations. In- 
deed, in the limit n — > oo, the effect of sea-quark pairs, 
such as the induced condensate or the OZI-rule violation 
in the vacuum channel, are enhanced. However, a large- 
Nf limit with N c held fixed would presumably meet phase 
transitions, leading eventually to the restoration of chiral 
symmetry and the loss of asymptotic freedom. The defini- 
tion of a large-fluctuation limit at the level of QCD is thus 
likely to require a combined limit of large Nf and large 
N c . However, at the level of the identities of the effective 
theory, we are allowed to take the formal limit n — > oo. 
We will see that this limit provides a natural example of 
large fluctuations, if we make the assumption that in this 
limit the pseudoscalar masses and decay constants remain 
finite order by order in powers of quark masses m and m s . 

Indeed, within such assumption, the fluctuation pa- 
rameters p and A grow with n, unless the constants L™ (/x) 
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and L 7 q((J.) are set to the asymptotic values of L^' "' and 



n,crit _ 



L 



where 



(80) 



Lf(/i) = 2Lf{») 



256tt 2 



M 2 
log— f- 



T log Mf 



( 81 ) 

From the chiral expansion of the masses and decay con- 
stants, Eqs. (|l7|)-(|lg|) and (p^)-(p^), and assuming that 
they remain finite order by order in powers of m and m s , 
one can read the a-priori unknown large-n behaviour of 
Y(n + 2), Z s n , £„, land A n , 



Y{n 



Zi 



£,n ~ l/l, A n ~ £ n ~ 1. 



(82) 



The constraint on Y(n + 2) is imposed by the presence of 
the logarithmic term (absent in the physical case n = 1): 



LL 



, M x 
log- 



-00 . 



(83) 



This implies in turn that Lg should grow like n, whilst L™ 
should approach a finite value, which need not coincide 
with the critical asymptote (pl|) . 

Actually, the same large-n behaviour of the theory is 
also imposed by the paramagnetic inequalities, provided 
that f n and g n , defined in Eqs. ( |42j ) and ( f45| ) are positive: 



fn ^ 9n 



M. 



2 r 



32tt 2 F 2 
nZoo > . 



, M x 
— log £ 

1 8 M 2 - 



log 



(84) 



The paramagnetic inequality for Z - second inequality 
in Eq. ( |63| ) - implies that y vanishes like l/n (or more 
quickly) . The limit y — > is allowed only on the first sheet 
and Eq. ( |60|) leads to v n — u n — ► oc like n 2 (or more 
quickly). On the other hand, the paramagnetic constraint 
for u - second inequality in Eq. ( |62| ) - leads tou = 0(n), 
unless a cancellation between kg and occurs. Such a 
cancellation between the NNLO remainder E and the one- 
loop chiral logarithms g cannot be logically ruled out. It 
would however stand against our initial assumption that 
the chiral expansions exhibit a good overall convergence 
starting at the NNLO order independently of n. 

Let us now describe more precisely the large-n be- 
haviour of u and v. Since u is positive, z = 1 — uy/A 
cannot diverge. Let us call Zoo = linin^oo z n and intro- 
duce a second parameter a such that: 



An, . 

(1 - Zoo) 

a 



4n 2 /a 2 , y~ a /n. (85) 



a and Zqc describe the behaviour of the fluctuation param- 
eters u and v. The order parameters tend to finite values: 



0. 



X(n + 2 
Z(n+2)^(l-ri(r)-e n )z 

X(2)[l-4„] 
Z{2)[l-e m ] 



1 - e(r) - 
[1 - n(r) 



(86) 



71 1 

e„][l 



akL 



In this limit, the LEC's associated with the decay constant 
identities remain scale-dependent - £5 = O(n) and L4 = 
0(1) - whereas the LEC's arising in the mass identities 
become: 



F 1 



IQrka 1 Ml ' 
' 1 — r\[r) 



(87) 



16a 2 



1 - e(r) - d n 



«r) + d' n ) 



Ml 



A comment is in order here about the double limit of large 
iV c and Nf, investigated in Ref. If we consider the 
large-n formulae for the LEC's IJ 56g [see Eqs. (|87|)-(j88|)l, 
we can recover their standard behaviour in the large- N c 
limit, provided that N c ~ n and a = 0(N C ), i.e. a = 
F 2 / A 2 with A an A^-independent scale 0. 

The large-n limit of the theory illustrates therefore 
how the large-fluctuation limit can be formulated consis- 
tently: the multi-flavour quark condensate is then sup- 
pressed whereas the two-flavour condensate A (2) remains 
different from zero (and close to 1 for r > 15). 



8.2 Large fluctuation parameters in three-flavour QCD 

The many-flavour limit of the theory discussed above should 
merely be viewed as a particular realisation of the limit 
of large fluctuations which is hopefully consistent but not 
necessarily unique. The limit of large fluctuations could 
as well be formulated directly in terms of fluctuation and 
order parameters keeping the number of flavours fixed to 
Nf = 3. In that case (n = 1), one avoids the presence of 
extra n 2 — 1 Goldstone boson X-states arising for n > 1. 
Since the limit is designed in a slightly different way from 
the previous section, we are ending up with similar but 
not identical results for the LEC's and the two- flavour 
parameters. 

One may infer from the generic Eqs. ( p4[ ) and ([35]) the 
lines in the (u, v) plane corresponding to a constant value 
of z = zi (0 < z\ < 1): 



y 



-(l-zr). 



r(i-y*i). 



(89) 



It is seen that keeping z fixed and setting y — > one repro- 
duces the large- fluctuation behaviour Eq. (|85|) . The decay 
constant Z(3) = [l—rj(r) — ei]zi remains nonzero, whereas 
the three-flavour condensate vanishes asymptotically: 



X(3) = [l- V (r)-e 1 }z 1 y^0. 



(90) 



In this large- fluctuation limit for n = 1, all relevant LEC's 
appearing in the original mass and decay constant iden- 
tities can be predicted in terms of r — m s /m and of the 



4 Needless to say that the "holomorphic phase" analysed in 
Ref. ]2£] is outside the scope of the present discussion which is 
exclusively based on an effective theory describing the break- 
down SU(Nf) x SU(N f ) -> SUv(Nf). 



S. Descotes, L. Girlanda and J. Stern: Chiral order and fluctuations in multi-flavour QCD 



13 



parameter Y(3) = X{3)/Z{3) 
of L$ and L% is then: 



0. The leading behaviour 



Y(3) ' 
e(r) + d[ 



8M 2 
Fl 



WM% Y 2 {3) 



(91) 



(92) 



Notice that the OZI-rule violating constant Le(fi) is no 
more suppressed [to be compared with the large-n result 
Eq. (p7|)]. This is best seen from the ratio L e /L s which in 
the y — > limit becomes: 



lim 



1 - e(r) 



(93) 



HisW (r + 2)[e(r)+4] ■ 

It amounts to: L 6 /L 8 = 0.43, 0.79, 1.53 for r = 20, 25, 30 
respectively. The second OZI-rule violating ratio L4/L5 
still depends on the order parameter Z(3): 



l im = 1 - ^ - Z ^ - ei 

v-oL B (ji) (r + 2)fo(r)+e£] 



(94) 



This ratio could be more suppressed provided Z{3) is close 
to 1. In the latter case the expression ([)i]) could be rather 
sensitive to the NNLO remainder ei, contrary to the case 
of the ratio Eq. ( p3| ) where d\ competes with 1 (recall that 
both d[ and e\ are suppressed by an extra power of 1/r). 

It remains to work out the two-flavour order parame- 
ters X(2) and Z(2) in the y — > 0, n = 1 large-fluctuation 
limit: 

X(2)[l-d\ a ] = (l-e(r)-d 1 )^—, (95) 



Z(2)[l-e„i] = (l-r?(r)- ei ) 



2 ' 
- 2zi 



(96) 



Despite the vanishing of the three-flavour condensate X(3) 
— > 0, the two-flavour condensate X(2) remains non zero 
and close to 1, provided r = m s /m is not too small. This 
effect is entirely due to the induced condensate and it 
is proportional to the strange quark mass. It is worth 
stressing the fundamental difference between the chiral 
symmetry restoration which occurs along the critical line 
v — M 2 /4 for finite u and v and the large-fluctuation limit 
in which (u,v) — > 00. This difference merely occurs for 
three-flavour order parameters: whereas in the symmetry 
restoration case both Z(3) and X(3) vanish and their ratio 
Y(3) = X(3)/Z(3) remains non-zero, the large fluctuation 
limit is characterised by a continuous decrease of Y(3) and 
of the condensate, with the decay constant Z(3) held fixed. 
This manifestation of large fluctuations need not corre- 
spond to a phase transition: they would lead naturally to 
a spontaneous breakdown of chiral symmetry with a very 
small (but non-vanishing) quark condensate. In this con- 
text, the recently proposed "no-go" theorems j2j| - stating 
that the vanishing of the condensate implies a vanishing 
pion decay constant and chiral symmetry restoration - do 
not necessarily apply. 

The limits of small and large fluctuations and of chi- 
ral symmetry restoration are summarised in Fig. ^, while 
the corresponding results for the order parameters are col- 
lected in Table liT 




Fig. 4. Limits of small and large fluctuations, and of chiral 
symmetry restoration on the y+ sheet for n = 1: the shaded 
region around the origin denoted N c corresponds to the domain 
of application for the large- N c limit, the arrow illustrates the 
large-fluctuation limit u, v — > 00, and the "x rest." border is 
the critical line of chiral symmetry restoration where z = 0. 



9 Two-flavour Chiral Perturbation Theory 

The striking outcome of the previous analysis is the per- 
sistence of a large two-flavour condensate S(2) even when 
the multi-flavour condensate H(Nf) [Nf > 3] is suppressed. 
We have seen that X(3) could be well below 1, indicating 
that the expansion of F%M% in powers of m u , rrid and m s 
need not be dominated by the genuine condensate term 
(m u + nid)S(3). Actually, the whole Nf > 3 xPT treated 
in the standard way can exhibit instabilities, not because 
of too large a strange-quark mass, but rather because of 
too sizeable vacuum fluctuations of qq pairs. On the other 
hand, the expansion of the same quantity F%M% in powers 
of m u , rrid only is expected to be dominated by the Nf = 2 
condensate term (rn u +md)S(2). Indeed, for not too small 
r = m s /m > 15, the result X(2) ~ 1 is emerging indepen- 
dently of the fluctuation parameters u and v reflecting an 
important contribution of the induced condensate. This 
suggests that for suitable pion observables, the standard 
SU(2) x SU{2) xPT g could be a well-convergent expan- 
sion scheme. 

In order to gain more insight into the different be- 
haviour of two- and multi-flavour chiral dynamics, it is 
convenient to rewrite the Nf — 2 Ward identities gen- 
erating the expansion of F%M% and of F% in a form as 
close as possible to Sees. || and f|. They involve the con- 
densate £(2), the decay constant F 2 (2) and the two 0(p A ) 
symmetry-breaking scale- independent LEC's £3 and I4 M: 



FlMl = 2mZ{2) 



m 2 B 2 (2) 
~ 8^ 



(4*4 



F 2 M 2 S (97) 
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Table 1. Values of chiral order parameters in three different limits (n = 1). The second column corresponds to small fluctuations 
(large- N c limit). The third one stands for the large-fluctuation limits (y — > 0, z set to z\). The fourth one indicates chiral 
symmetry restoration (z — ► 0, y set to j/i). 



Limit 


AT , 
iV c — > 00 


Large fluct. (y — > 0) 


X rest. (2 — > 0) 


Ml 

vi 






4(1 -yzi)/?/ 2 


4/2/1 

4/y? 


Y(3) 


1 - e(r) - di 
1 - r]{r) - ei 





1 - e(r) - di 

1 n ^ 

1 - i](r) - ei 


Z(3) 
Z(2) 


1 - r]{r) - ei 
1 — r](r) - ei 


zi[l — r?(r) - ei] 
r + 2jzi 1 — ??(r) — ei 




r 1 - 77 (r) - ei 


r + 2 1 - 


r + 2 1-^ [1 ^ lfcl(r)3l] 


X(3) 
X{2) 


1 - e(r) - di 
1 - e(r) - dj 




r 1 — e(r) — di 




r 1 - e(r) - di 2; , w 1 


r + 2 1 - eki 


T + 2 1 — d,rl 



F 



F 2 (2) 



mB(2) 
4tt 2 



F 2 £ . 



(98) 



Here S(2) = Z"(2)/F 2 (2) and the NNLO remainders S 
and e [] are 0{m 2 ), expected to be of order 1%. The anal- 
ogy with the multi-flavour case can be pushed further by 
rewriting Eqs. @ and @ in the form of Eqs. ©-©: 



Z{2) = 1 - e- Y(2)A/4, 
(99) 



X(2) = l-5-Y(2) 2 p/4, 
where: 

1 M 2 - 
The simple rescaling: 



1 M 2 - 

(100) 



1 F 2 (2) 



y 



l-s 2mB(2) 



1-e F 2 ' a l-<5 
u = Xk , v — pk , k — 



Ml 

1-6 
(1-e) 2 



(101) 



(102) 



brings the two-flavour mass and decay constant identities 
into the standard form ( 54 )- (p5j) : 



z + -uy 



-,vy 



1, 
1 

y 



(103) 
(104) 



The formal analogy is now complete. In the two-flavour 
case the factors e(r) and 77(f) have to be omitted in the 
rescaling factors but otherwise, all equations look identi- 
cal. Where is the difference? 

The multi- and two-flavour cases behave differently be- 
cause the corresponding parameters v and u (or p and A) 
are respectively of a different origin and magnitude. We 
know from the previous discussion that the two-flavour 
quantity y ~ 1 and consequently, the parameters v — u and 

5 The remainder e should not be confused with the function 
e(r) introduced in Section 3.2. 



u cannot be too large. This in turn excludes an unlimited 
grow of |^3 1 . As a consequence one may expect the per- 
turbative solution of the fundamental equations (103) and 
(104) [obtained by Taylor expanding the non-perturbative 
solution] to make sense. 

We have already emphasised the connection of the 
multi-flavour parameters v n ,u n (n > 1) with the corre- 
lations between vacuum strange and non-strange qq pairs 
and (last but not least) with the fluctuations of small Eu- 
clidean Dirac eigenvalues H . We do not know much about 
the importance of these fluctuations from first principles, 
but we do understand why in the Nf > 3 theory such fluc- 
tuations manifest themselves through important OZI-rule 
violations in the vacuum channel J p — + . Hence, if the 
quark pairs in the vacuum are strongly correlated and/or 
the low-energy Dirac spectrum is subjected to large fluctu- 
ations, the multi-flavour parameters u n , v n are likely large 
and the perturbative solution of the corresponding fun- 
damental equations ( |54| ) and (|||) breaks down. On the 
other hand, in a Nf = 2 theory and in the presence of 
massive (m s ~ Aqcd) strange quarks in the sea, the same 
cause does not produce the same effect. In this case, the 
fluctuations of small Dirac eigenvalues are much harder 
to relate to low-energy observables: the OZI-rule is inop- 
erative in this case, and the scalar correlator ({uu)(dd)) 
is chirally invariant and not simply related to an observ- 
able order parameter. The different nature of parameters 
u and v is further illustrated by a different behaviour in 
the large- N c limit: whereas the multi-flavour fluctuation 
parameters u n , v n for n > 1 are suppressed as 0(1/N C ), u 
and v behave as 0(1) since the constants £.3 and £4 behave 
like 0(N C ). 

The fact that v, u (or p, A) as well as the LEC's £3,^4 
need not be enhanced by fluctuations of small Dirac eigen- 
values can be seen directly by comparing (matching) the 
multiflavour mass identity (|17| ) with its Nf — 2 counter- 
part (pT^). In the former, the vacuum-fluctuation contri- 
bution resides in the term containing the LEC Lq(p). 
Eq. d97| ) is a reexpression of the same mass identity in 
terms of 5(7(2) x 5(7(2) quantities. Both identities have 
to match order by order in m. At the first order, one ob- 
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tains the expression for the two-flavour condensate: 



17(2) = £(n + 2) 



(105) 



which absorbs the major part (the term proportional to 
m s ) of the vacuum fluctuation contribution (nm s +2m)Z s . 
nm s Z^ is what we have called the induced condensate j|, 
|24|. The remaining part of Z s contributes together with 
A n (i.e. L s ) to the Nf = 2 LEC's £3 and £4, but this contri- 
bution contains neither m s nor the flavour factor n. Con- 
sequently, there is no particular reason why fluctuations 
of small Dirac eigenvalues should enhance the parameters 
v and u and to destabilise the two-flavour standard xPT. 

Th e no n-perturbative effects in the non-linear system 
( |l 3| ) - ( |T04j ) start to show up for \v — u\ ~ 1 correspond- 
ing to |^3 1 ~ 35. At this point, it is worth stressing that 
the original standard estimates of £3 , £4 make use of 

Nf — 3 observables assuming the approximate validity of 
the OZI-rule in the scalar channel. Since vacuum fluctu- 
ations do contribute to £3, I&, the standard estimates ||, 
0] £3 = 2.9 ± 2.4 and £4 = 4.3 ± 0.9 could be modified by 
OZI-rule violating effects. 

In order to relate £3, £4 to Nf = 3 observables includ- 
ing OZI-rule violating vacuum fluctuations, one may pro- 
ceed as follows. First the two-flavour mass and decay con- 
stant identities Eqs. ( p7| ) and (|9^) are equivalently rewrit- 
ten to express £3, £4 in terms of the order parameters X(2) 
and Z{2) : 



£3 = 32tt 



2 F*Z{2) 
M%X{2) 



1-e - (1-5) 



Z{2) 
X{2) 



, (106) 
(107) 



Next, one uses Eqs. @ and @ relating X(2) and Z(2) 
to the multi-flavour order parameters y n and z n . In the 
case of three flavours (n — 1) the latter may be written 
as: 



Z{2) 



X(2) 



r + 2 



[1 - r)(r) - (e x - e ni . . .)] 



x \ 1 - yih(r)gi + -zi> , 

[1 - e(r) - (di - d wl . . .)] 



r + 2 



x { 1- yfk 1 (r)f 1 + -y t zi 



(108) 



(109) 



where the ellipses stand for terms of order (e^i) 2 , (d^i) 2 
as well as e^i^r) and d„ie(r). Since e^i, d n i = O(rrij ) are 
expe cted of order 0.1 or less, the dotted terms in Eqs. (108) 
and (|109| ) can hardly exceed 1 %. Furthermore, the leading 
0(mj.)_ terms cancel out in the differences e\ — e„i and 
di — d^i- As a resu lt, t he wh ole c ontribution of NNLO 
remainders in Eqs. ( |108] ) and (109) is 0(rnm s ), i.e. they 
can be expected of order 3 %. 

It is instructive to exploit the above relations and to 
estimate the order parameters X(2),Z(2) and the LEC's 



£3, £4 that correspond to the extreme cases of no fluctua- 
tion (N c — > 00) and large fluctuations (yi — > 0). Neglect- 
ing all NNLO remainders, one gets for N c 00: 



20 



r = 25 



X{2) = 0.905 
h = -7.5 , 



Z(2) = 0.949, 
£4 = 2.0 , 



X(2) = 0.955, Z(2) = 0.959, 
£3 = -0.6 , £4 = 1.5 , 



(110) 



(111) 



These estimates are quite compatible with standard ex- 
pectations [BIQ]. On the other hand, one gets in the limit 
of large fluctuations (taking Z\ = 0.7): 



20 



r = 25 



X{2) = 0.822 , Z{2) = 0.922 , 
£3 = -20.2, £4 = 3.2, 

X{2) = 0.884, Z(2) = 0.938, 
£3 = -9.6 , £4 = 2.4 , 



(112) 



(113) 



As expected, large vacuum fluctuations push £3 to- 
wards larger negative values, without however reaching 
the range in which the two-flavour xPT would start to di- 
verge. The above estimates should be ta ken w ith c aution: 
even the small remainders in Eqs. ( |l08[) and (|lQS| ) of the 
order of a few per cent can give an important contribu- 
tion to £3 and £4 since the latter involve X(2) — Z(2) and 
1 — Z(2) respectively. 

The information on low-energy nn phases extracted 
from the new E865 Brookhaven K e 4 experiment | ft5| com- 
bined with older data on I = and 1 = 2 S-waves p6[ , 
have been recently used to extract the two-flavour order 
parameters X(2),Z(2) as well as the LEC's £3,^4 in a 
model-independent analysis [po| which is merely based on 
the numerical solution of Roy equations Jl^|. The result 
of this analysis reads: 



X{2) = 0.81 ±0.07, 
£3 = -17.8 ± 15.3, 



Z{2) = 0.89 ±0.02, (114) 
£4 = 4.1 ±0.9, (115) 



It is compatible with the estimates given above, and espe- 
cially with the ones corresponding to the large-fluctuation 
limit. 

We would like to close this section with a comment em- 
phasising the importance of the new precise 7T7r-scattering 
data recently published Jl5| , forthcoming (2?J or expected 
in a near future J28|. Some time ago it has been pointed 
out that no experimental test of the actual size of the 
quark condensate and of the magnitude of the quark mass 
ratio r = m s /m was available. No particular mechanism 
of xSB was anticipated at that time and a general (less 
predictive) expansion scheme (Generalized Chiral Pertur- 
bation Theory or G%PT) was proposed to analyse data in 
a model-independent way. Furthermore, there was no or 
little experimental evidence for the substantial violation 
of the OZI rule in the scalar channel that supports today 
the idea of important vacuum fluctuations and a sizeable 
difference between the two- and multi-flavour condensates. 

On the other hand, the present conclusion that X(2) 
is close to 1 and its consequences for the standard two- 
flavour X-PT might appear at first sight as a result that 
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we could have anticipated using nothing but theoretical 
arguments similar to those in this paper - contradicting 
therefore the claims made in Ref. [£9|. As a matter of 
fact, this conclusion can be drawn only if one gets ex- 
tra (and independent) information about the size of the 
quark mass ratio r = m s /m, excluding small values such 
as r ~ 10. In the present paper, a larger value of r was as- 
sumed from the onset. If r happened to be small, the factor 
e(r) = 2(r2 — r)/ (r 2 — 1) would get close to 1, affecting the 
resc aling factors in Eqs. (|52[)-(p3|) and (|56|). According to 
Eq. (10E), X(2) would then be suppressed independently 
of the size of vacuum fluctuations. Because of the inequal- 
ities X{2) > X{3) > 0, both X{3) and X(2)-X{3) would 
then be suppressed. Hence, for small r, the vacuum fluc- 
tuations and the induced condensate would enhance X(2) 
less significantly, implying a weaker violation of the OZI 
rule in the scalar channel. In all cases, a strong, correlation 
between r and X(2) [not X(3)] persists [Q,|l3|]. Hence an 
important result of the new and/or forthcoming 7T7r scat- 
tering experiments is - among other issues - to put a lower 
bound on the quark mass ratio r, ruling out small values 
of the two- flavour GOR ratio X{2). A precise quantita- 
tive statement of this lower bound is matter of a detailed 
SU(3) x SU(3) analysis of 7T7t scattering along the lines of 
the present paper. Such an analysis is in progress 130]. 



10 Conclusion 

i) Chiral order parameters may exhibit a strong de- 
pendence on the number of light flavours. Such a phe- 
nomenon can be interpreted as an important effect of sea- 
quark pairs, which is related to a large violation of the 
OZI-rule observed in the scalar channel and reflects sig- 
nificant fluctuations of the lowest modes of the Euclidean 
Dirac operator. In this paper we have focused on the pre- 
cise interplay between chiral order and fluctuations. 

ii) We have highlighted the particular role of the strange 
quark, whose mass m s ~ Aqcd is light enough to popu- 
late the vacuum with massive ss-pairs, but heavy enough 
to have influence on the properties of the SU{2) x SU(2) 
chiral limit with m u , md —* and m s fixed at its physical 
value. In particular, the two-flavour condensate S(2) = 



lim„ 



,^,o( uu ) does receive an extra contribution from 



the massive vacuum ss-pairs through the OZI-rule vio- 
lating and SU (2) x SU (2) symmetry-breaking correlation 
((uu)(ss)). The latter is referred to as the induced conden- 
sate. Such contribution is absent in the SU{3) x SU(3) 
chiral limit m u ^d, s — * 0, since the remaining quarks (c, 
b, t) are too heavy to contribute significantly to vacuum 
fluctuations. This may lead to significant differences in the 
manifestation of chiral symmetry breaking in the two- and 
three-flavour chiral limits. 

Hi) To probe the effect of massive quark pairs on chi- 
ral symmetry breaking, we have introduced a theory with 
two ultralight u, d quarks and n degenerate copies of the 
(massive) strange quark. QCD can be recovered as the 
particular case n = 1. Because of the mass hierarchy 
mu,d <C m s <C Ah, two different chiral limits can be con- 
sidered in this model: the two-flavour limit, where only 



niu,d — > but m s ^ 0, and the multi- flavour limit, where 
all Nf — n + 2 light masses vanish. 

iv) The multi- flavour case is characterised by two fluc- 
tuation parameters which measure the violation of the 
OZI-rule in the scalar sector. They are related to the two 
large- N c suppressed 0(p 4 ) constants L 6 ([i) and L±(n). We 
have shown that Ward identities yield two non-linear rela- 
tions between the fluctuation parameters and the two ba- 
sic multi-flavour order parameters: the quark condensate 
X(N f ) = 2mU(N f )/(F^M^) and the Goldstone-boson 
decay constant Z(N f ) = F 2 (N f )/F% (where F(N f ) is the 
pion coupling constant F% in the chiral limit). These rela- 
tions are a direct consequence of Ward identities: all higher 
chiral orders (NNLO and beyond) are absorbed into a fi- 
nite multiplicative renormalisation (rescaling) of the order 
and fluctuation parameters. The effect of this renormali- 
sation remains small (i.e. rescaling factors ~ I) provided 
that the chiral series globally converge from the NNLO 
order. 

v) Taking the SU(2) x SU(2) limit m u ,m d -> of 
mass and decay constant Ward identities, one obtains the 
two-flavour condensate X{2) and decay constant Z(2) in 
the presence of n copies of massive s-quark pairs in the 
sea. We can then compare multi- and two-flavour order 
parameters as functions of fluctuation parameters: 

— Multi-flavour chiral symmetry is restored along a crit- 
ical line in the plane of fluctuation parameters. Along 
this line, the multi-flavour condensate X{N f) and de- 
cay constant Z(Nf) both vanish, but their ratio X/Z 
stays finite and non-zero. The two-flavour condensate 
and pseudoscalar decay constant do not vanish in this 
limit, except in one exceptional point corresponding to 
the end-point of the critical line. 

— In the case of small fluctuations, we recover a large- N c , 
mean-field type behaviour. The order parameters then 
do not depend on the number of light flavours, and 
the quark condensate is the dominant signal of chiral 
symmetry breaking. 

— In the opposite limit of large fluctuations, the multi- 
flavour quark condensate X(Nj) tends to zero, but 
chiral symmetry remains spontaneously broken, since 
the decay constant Z(Nf) stays away from 0. Corre- 
spondingly, the two-flavour order parameters receive 
large contributions from massive sea-quark loops, so 
that their size is naturally large, and close to the one 
expected in the large- N c limit. We would like to stress 
that the presence of large fluctuations is not neces- 
sarily related to the occurrence of a phase transition: 
they could constitute a feature of the dynamics of the 
theory, implying in turn a naturally small size of the 
quark condensate. In this case the two-flavour conden- 
sate would be exclusively made of the induced contri- 
bution from the massive quark pairs in the vacuum. 

vi) The usual treatement of Nf > 3 x?T considers 
fluctuation parameters as small, which requires a very pre- 
cise fine-tuning of the 0(p 4 ) constants Lq and L4. On the 
other hand, the available information about the OZI-rule 
violation in the scalar channel suggests that fluctuations 
can be enhanced due to the large effect of sea quarks. 
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Such large fluctuations destabilise chiral expansions by 
suppressing the lowest order of the series (three-flavour 
quark condensate), but they need not spoil the overall 
convergence. To cope with such a situation, we propose 
to treat the fluctuation parameters (related to the OZI- 
rulc violating LEC's L4 and Lq) non-perturbatively when 
expressing the parameters of the effective Lagrangian in 
terms of observables. In practice it amounts to replac- 
ing the perturbative solution of the non-linear relation 
between order and fluctuation parameters (expressed as 
an expansion in powers of the latter) by the correspond- 
ing exact algebraic solution. In this way one can reex- 
press the low-energy constants tuBq, F , L 5 , L7 and L 8 in 
terms of the two fluctuation parameters, the quark mass 
ratio r — m s /m and higher xPT contributions arising 
through the (hopefully convergent) expansion of rescaling 
factors around 1. This procedure corresponds to a non- 
perturbative resummation of the fluctuations encoded in 
L4 and Lq. 

vii) We have shown that vacuum fluctuations do not 
suppress the two-flavour condensate. The Nf = 2 GOR ra- 
tio X(2) remains close to 1 and the standard two-flavour 
xPT is likely a well-convergent expansion scheme, pro- 
vided that the quark mass ratio r is not too small. A 
lower bound on r can be inferred from precise low-energy 
7T7r scattering data. 

viii) Additional observables (Goldstone boson scatter- 
ing, decay form factors. . . ) would have to be analysed 
along similar lines, considering their chiral expansion up 
to NNLO and making explicit the non-perturbative de- 
pendence on the fluctuation parameters. We will illustrate 
more precisely this point for tttt and irK scattering param- 
eters in a future publication Q. Applying this method 
to a sufficiently large set of precisely measured observ- 
ables should allow one to pin down the size of vacuum 
fluctuations, to disentangle the effect of massive ss-pairs 
on the pattern of two-flavour chiral symmetry breaking 
and to determine by how much three-flavour chiral order 
parameters are suppressed. This will eventually lead to 
a better understanding of the spontaneous breakdown of 
chiral symmetry and its dependence on the number and 
hierarchy of light quarks. 
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acknowledges partial support from European Program HPRN- 
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A Goldstone boson masses in multi-flavour 
QCD 

We want to discuss in this section how to combine the 
Ward identities for the masses and decay constants of the 
Goldstone bosons. We are going to see that we need as 



inputs M w , Mk, F n , Fk and M v to fix the values of the 
0(p 4 ) LEC's as functions of r, B 0n and F 2 (n + 2), and 
then to obtain the values of the remaining observables F 2 , 
Mf and F x . 

First, we have to spell out all the Ward identities con- 
cerning the masses and decay constants of the Goldstone 
bosons. The masses and decay constant identities for the 
pion and kaon were given in Eqs. ([l7])-([l8]) and Eqs. (p5|)- 
(p(f). The Ward identities for the extra X states can be 
expressed in the more practical form: 

FlMl 



F 2 M X 
FlMl 



4r 



hi 
4r 



F 2 K M 2 K 



(116) 
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F 2 K 
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F 2 M 2 
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F X M\ 



Xn , 



16tt 2 F, 



A 



1 M 2 
2 l0g Mf 



r M 2 
-log— £ 

n g M 2 X 



F 2 k 



exn 



-2e K 

■'I, 

F K 



(117) 



We have not discussed yet the identities for the 77, which 
can be recasted in a form reminiscent of the Gell-Mann- 
Okubo formula: 

(n + 2)F 2 M 2 - AF 2 K M 2 K - (n - 2)F 2 M 2 

= 4(r-l)m 2 {(r-l)(2nZP + A n ) 

+B 2 n [(r - l)L' n - L}} + (n + 2)F 2 M 2 d v 
-AF 2 K M 2 K d K - (n - 2)F 2 M 2 d 7T , (118) 



(n + 2)F 2 - 
2m_Bor, 



AF 2 K 



16tt 2 



(n - 2)F 2 

— log —J- - (2n - 1) log —I 



M 2 X 



-2(7^ -1) log -f 



{n + 2)F 2 e v - AF 2 e K 



-{n-2)F 2 e n . 



(119) 



The ry-mass identity involves the new LEC = \%B\ n LVj . 
We have also introduced the NNLO remainders d v and e v . 

The Ward identities are therefore expressions of Fp 
and FpMp (P — tt, K, rj, X) in terms of the fundamen- 
tal parameters r = m s /m, S(n + 2), F 2 (n + 2), the LEC's 
L4...8, chiral logarithms of pseudoscalar masses and NNLO 
remainders. The pion and kaon identities yield then £4,5,6, 8 
in terms of r, Bo n , F 2 (n + 2), chiral logarithms and NNLO 
remainders, see Eqs. @), (§|), @ and @. We can 
obtain a similar expression for Lj from the identity for 
F 2 M 2 . From Sec. [|, pion and kaon Ward identities lead 
to the expression of the parameters -Bon and F 2 (n + 2) as 
functions of the two fluctuations parameters u, v, the ratio 
of quark masses r and NNLO remainders [see Eqs. (|54|)- 
©]■ 

We can now exploit all the remaining Ward identities 
to write F x , M x and F 2 as functions of r, u, v, chiral log- 
arithms of Mp and NNLO remainders. We need thus to 
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Table 2, 

n = 2. 



Value of Mx/Mk as a function of Y(4) and r for 



Y(4) 


r = 20 


r = 25 


r = 30 





1.43 


1.37 


1.30 


0.2 


1.44 


1.38 


1.31 


0.5 


1.44 


1.38 


1.31 


1 


1.46 


1.40 


1.32 


2 


1.48 


1.42 


1.35 


4 


1.57 


1.53 


1.46 



Table 3. Mass logarithms defined in Eq. (121) and involved in 
Z{2) and X(2), for r = 25, n = 1 and Z(3) = 1 (j = 0.021 and 
hi = 0.059). M nt K,r] and F^x are set to their physical values, 
and all NNLO remainders are neglected. 



Y 


6 


c 


/ 


9 





-0.018 


-0.021 


0.037 


0.062 


0.5 


-0.013 


-0.015 


0.043 


0.067 


1 


-0.008 


-0.009 


0.050 


0.072 


2 


-0.003 


-0.002 


0.056 


0.077 


4 


-0.009 


-0.005 


0.053 


0.071 



know the masses M Wl K,ni the decay constants F w ^k and 
the NNLO remainders. As an illustration of the method, 
we take the physical values of the masses and decay con- 
stants and set the remainders to zero, in order to obtain 
M x through an iterative procedure. We start with the 
approximate value: 



Ml 



rF 2 



2F 2 

K 



F 2 



FlMl 



1 



(120) 



and iterate Eqs. ( |116| )- ( J117| ) until they converge to M x - 
For instance, the values for Mx /Mk a.tn — 2 are collected 
in Table @. 



Tab le 4 . Combination of chiral logarithms defined in 
Eq. dl27 ) and involved in Z(2) and X{2) when n —> oo. M^ yK 
and i^if are set to their physical values, and all NNLO re- 
mainders are neglected. 



a 


r = 20 


r = 25 


r = 30 





0.003 


0.003 


0.002 


1 


0.003 


0.003 


0.002 


2 


0.003 


0.003 


0.002 


4 


0.004 


0.004 


0.003 


6 


0.004 


0.004 


0.004 


10 


0.006 


0.007 


0.007 



B Masses in the ST/ (2) x SU(2) chiral limit 

To know the masses of the Goldstone bosons in the SU (2) x 
SU (2) chiral limit, we take the limit m — > in all the mass 
and decay constant identities, and reexpress each LEC in 
this limit in terms of the same LEC with and chiral 

logarithms of Mp/Mp, see E_qs. @ and Q. Using the 
previous relations, Fp and Mp are functions of (r, u, v), 
pseudoscalar masses Mq and Mq, and remainders. If we 
keep on setting M^^.-q and F^^r to their physical values, 
and neglecting NNLO remainders, we can compute all the 
masses in the SU(2) x SU(2) limit in terms of the ratio of 
quark masses r and the fluctuation parameters u and v. 

We need the pseudoscalar masses Mp in particular to 
compute the factors /„ and g n arising in the expression 
of SU(2) x SU(2) order parameters Eqs. @ and @. 
These factors can be written as: 



Ml r , 



Mlr+1 1 



fn 



n 2 



M 2 
i v 

JJT^l^ 2 ° g M|' 

M 2 

= --f(nr + 2)^„, 
h n + j + c n . 



(121) 

Table H shows an illustration of the physical case n = 1 . 
The quark mass ratio is set to r = 25. Z(3) and Y(3) were 
chosen as parameters, rather than uj and v\ [the two sets 
are related through Eqs. ( pi|) and (p5|)], and we set Z(3) 
to 1. j is only a function of r. h n depends on Y(n + 2) only 
through L' n , which vanishes for n = 1. For r = 25, we have 
j — 0.021 and hi = 0.059. Furher study shows that fx and 
gi are only weakly dependent on Z and r. If we set the 
masses and decay constants to their physical values and 



vary n, we observe that /„ 
of n. 



and g n are increasing functions 



C Large-n limit 

We take now the large-n limit as described in Sec. ^.1[ 
We stress that this limit is considered only at the level 
of the effective theory, with the additional assumption 
that the masses and decay constants tend to some par- 
ticular asymptotic values. For numerical purposes, we will 
furthermore take these asymptotic values as the physical 
ones, and set the NNLO remainders to 0. We recall that 
we have introduced in this limit two parameters and a 
describing the behaviour of SU (n + 2) x SU(n + 2) chiral 
order parameters: 



Z(n + 2)^(l-n(r)-e n )z 
Y(n 



2 - al- e(r) - d n 



(122) 



n 1 — r)(r) 



The chiral logarithms disappear then from the mass iden- 
tities Eq s. (|l^)-([ig|). In the decay constant identities Eqs. 

only the logarithms of M 2 K jM 2 x survive - M v 
disappears from the right hand side of all mass and decay 
constant identities. Due to this simplification, Eq. (120) 
becomes an exact formula for Mx (up to NNLO remain- 
ders). For r = 20,25,30, M X /M K is respectively 1.43, 
1.37 and 1.30. 

As far as the 77-meson is concerned, the mass iden- 
tity Eq. (118) and the finiteness of all masses and de- 
cay constants at large n result only in a constraint on 
the large-n behaviour of Z? (i.e. L"). On the other hand , 
the (LEC-free) relation for the decay constants Eq. ( 119 ) 
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imposes that — F v up to NNLO remainders when 
n — » oo. This can be related to the structure of the rj 
meson: A r) = \J\ + 2/n- diag[l, 1, —2/n. . . — 2/n], so that 
a similar equality may apply to tt and 7/ masses in this 
limit. Such an assumption is however not necessary for 
our purposes: as we noticed earlier, the large-n behaviour 
of M v does not affect the K and X-meson spectrum since 
\ogM 2 disappears from the Ward identities for the un- 
mixed states. 

Following the same lines as in the previous section, we 
can exploit the Ward identities in the SU(2) x SU (2) chiral 
limit to determine the pseudoscalar masses for m — ► 0. Let 
us in particular notice the identities for the kaons: 
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In Eq. (12E), we would naively expect the first line to 
be dominant at large r. This is not actually the case. We 
have: 



r F 2 M 2 r 



2F 2 K M 2 K 



+ 1 ' 



r% = 2 



M 



K 



Ml 



1 = 24, 



(130) 

and (e — rf) (r) is precisely changing sign at , being pos- 
itive for r < r\ and negative r > r\. We see ther efore 
that the two logarithms in the first line of Eq. (|129| ) are 
of opposite signs, exchanging them for r ~ 24. 

The actual sign of is therefore a question of subtle 
compensation between all the logarithms involved in its 
expression. Following the procedure outlined above, we 
have computed the values of for various values of r 
and a, collected in Table ||, with M ff) x and F^.k set to 
their physical values, and all NNLO remainders neglected. 
We see that remains small and positive in any case. 
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Comp_aring Eqs. (|125|)-(|l26|) with Eqs. (|ll6[)-(|ll7|) shows 
that F x — F x and M x = M\ in the large-n limit (up to 
NNLO remainders). 

The chiral logarithms involved in the discussion of X(2) 
and Z(2) reduce now to a single combination of loga- 
rithms, independent of rj observables: 
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If we neglect (as a first approximation) the chiral log- 
arithms in Eq. (122) - setting a = - we obtain: 
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